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Abstract 

We show that geometric confinement dramatically affects the shear-induced configurations of 
dense mono-disperse colloidal suspensions; a new structure emerges, where layers of particles buckle 
to stack in a more efficient packing. The volume fraction in the shear zone is controlled by a 
balance between the viscous stresses and the osmotic pressure of a contacting reservoir of unsheared 
particles. We present a model that accounts for our observations and helps elucidate the complex 
interplay between particle packing and shear stress for confined suspensions. 
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Colloidal suspensions in thermodynamic equilibrium exhibit fascinating phase behavior, 
controlled by a delicate interplay between inter-particle interactions and volume exclusion. 
Packing constraints are essential for determining the structures that result and an under- 
standing of these has been instrumental in elucidating the phase behavior formed by colloidal 
suspensions under various conditions. These include the packing of spherical particles in bulk 
Q, 0| , and the unavoidable defects formed when particles order on the surface of a spheri- 
cal drop 0| . Exposing such suspensions to large strains can drive them out of equilibrium 
and significantly modify the particle configurations; nevertheless, the shear stresses can still 
effectively thermalize the particles, allowing them to explore their phase space, and adopt 
reproducible structures. For example, in bulk, a dense suspension of mono-disperse particles 
subjected to oscillatory shear wi. 
parallel to the shearing plates 
layers to separate allowing them to flow over one another, and the particle velocity and the 
amplitude of the motion vary linearly between the shearing plates . While these 

structures have been well described, the complex interplay between particle packing and the 
shear induced stresses which leads to formation of these structures is still poorly understood. 

When a dense suspension of mono-disperse particles is geometrically confined between 
two plates but not subjected to shear, the packing constraints force the suspension to adopt 
equilibrium structures different from those observed in bulk 



l_order into hexagonally-close-packed (HCP) layers oriented 
E|. The shear-induced viscous stresses force adjacent 
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1 it ] . Subjecting such 



highly confined suspensions to shear is of considerable technological relevance to coatings, 



lubricants, and bio-rheology 
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141 ]; moreover, the limited number of particle layers may 



facilitate a more quantitative analysis of the stresses, and a direct determination of the in- 
terplay between packing and shear-induced stresses, allowing the resultant non-equilibrium 
structures to be explicitly determined. Surprisingly, such highly confined suspensions sub- 
jected to shear have never been investigated. 

In this paper, we investigate dense colloidal suspensions highly confined between two flat 
plates and subjected to large oscillatory shear. We show that confinement forces the sus- 
pension to adopt structures that are translationally invariant along the direction of particle 
motion, and include striking gaps in the packing which, nevertheless, allow the particles to 
pack more efficiently than those observed in bulk. We present a model that accounts for 
our observations by elucidating the interplay between shear stress, particle packing, and 
geometric confinement that leads to these ordered, but highly non-equilibrium structures. 
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FIG. 1: Structure of a sheared suspension with 7 = 0.38, / = 30Hz, and = 0.61. The plate 
moves in the x direction. Figure la shows a confocal micrograph of a sheared suspension forming 
HCP layers when D = 80fim. Figures lb-e show micrographs of the suspension in the buckled 
state. The gap is set slightly below the height commensurate with confinement of four flat HCP 
layers. Figure lb shows an x-y image slice of the suspension near the upper plate. Figures lc, 
Id, and le show slices that are, respectively, 1.3/iin, 2.6/mi, and 3.9/um below the slice in lb. The 
images are presented side by side so that the strip alignment can be compared. The y' direction is 
aligned with one of the characteristic HCP lattice vectors and forms a 60 degree angle with the x 
direction. The x, y, and y' directions are indicated in the bottom left corner of la. 

The suspensions contain spherical poly-(methylmethacrylate) particles sterically stabi- 
lized by a thin layer of poly-12-hydroxystearic acid The particles are impreg- 
nated with rhodamine dye and suspended in a mixture of cyclohexyl bromide, tetralin, and 
decalin, chosen to match both the index of refraction and density of the particles. This 
procedure allows us to view the three dimensional structure of the suspension using fluores- 
cence and confocal microscopy 16]. Optical tweezers measurements |l8j] and electrophoretic 
mobility measurements show that the dying process imparts a positive charge onto the par- 
ticles. The liquid- crystal coexistence regime for these particles is shifted and has been shown 
to occur at a particle volume fraction 0.38 < <p < 0.42 [19] whereas in hard spheres, this 
regime occurs at 0.494 < <fi < 0.545. While the experiments described in this paper use 
charged particles, preliminary experiments using particles that are significantly less charged 
reproduce the observed phenomena and indicate that charging effects play a secondary role 
in the observed pattern formation. The solvent mixture viscosity, rjo = 0.023poise. The 
particles have a diameter d = 1.42/im and polydispersity of 5%. 

Our apparatus jscj mounts on an inverted microscope and allows 3-D visualization of 
a suspension with control over gap separation, shear frequency, and shear amplitude. The 
suspension is sheared between a fixed plate 5mm in diameter and a movable microscope 
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cover slip, both of which are flat on the particle length-scale. Set screws fix the gap, D, 
between 1/mi and 100/zm, and align the plates so they are parallel to within one micron 
over the shear zone. A sinusoidally driven piezoelectric actuator moves the cover slip and 
produces up to 30/zm displacements at frequencies, /, ranging between 0Hz and 100Hz. 
The 5mm plate is immersed in the suspension so particles in the shear zone contact a large 
reservoir of particles fixed at volume fraction, = 0.61 ± 0.02. The cell is enclosed so no 
solvent evaporates. This lets us work with the suspensions for periods longer than a year. 

To drive the suspension out of equilibrium, we shear with maximum strains, 7 > 0.3 and 
frequencies, / > 5Hz. A phase diagram of the transition from equilibrium to nonequilibrium 
structures as a function of 7 and / will be shown elsewhere 211]. When the gap holds 
more than 11 layers, the morphology is identical to that in bulk suspensions (Fig. |TJi). 
However, when the gap holds fewer than 11 layers, confinement plays a critical role. The 
HCP layer structure becomes intermittent, occurring only at discrete plate separations. 
For gaps incommensurate with these separations, we observe a remarkable new ordering 
(Figs[2>e). Fluid voids appear within the planes and the layers break up into HCP strips 
aligned in the direction of plate motion, x. Moreover, the strip widths vary with depth, z. 
A typical example obtained for a gap slightly smaller than that confining four flat layers, 
is shown by the series of x-y images in Figs [TJ>e. Near the stationary upper plate, the 
strips are three particles wide and have fluid voids between them (Fig. Gfc). In a plane 
1.3/im lower, particles orient in two-particle-wide strips, alternating with one-particle-wide 
strips (Fig. Efc). Remarkably, the velocity of the two-particle-wide strips is larger than 
that of the one-particle- wide strips. This structure repeats 1.3/zm further down, but this 
time the velocity of the one-particle-wide strips is larger (Fig. [TJl). Finally, 1.3/im further 
down, the layer nearest the bottom shearing plate is again oriented into three-particle-wide 
strips separated by fluid voids (Fig. [T^). It is convenient to examine the strips along the 
HCP lattice vector directions. Inspection of Fig. 0d along the y' axis shows that the three- 
particle-wide strips are registered. Microscopy measurements indicate the registration arises 
from interdigitation with the one-particle-wide strips which align with the fluid voids but 
are located in the plane 1.3/im below (Fig. [It). This configuration forces the one and three 
particle-wide strips to have equal velocities. Similar interdigitation is observed for the two- 
particle- wide strips in Figs^ andCJi. Again, the interdigitation forces the strips in different 
layers to have equal velocities. Finally, similar behavior is observed for the one and three 
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particle- wide strips in Figs[T]i and[T^. 

We summarize this behavior in the y'-z schematic of the particle positions shown in 
Fig. |21 The dashed lines indicate the z position of the x-y planes depicted in Figs[T)>e. The 
shadings delineate particles in contact and moving with equal velocities. This figure shows 
the peculiar patterns result from a buckling of the colloidal and lubricating fluid layers. The 
particle velocity and oscillation amplitude of the buckled particle layers still vary linearly 
between the shearing plates. Since the two sets of particle strips shown in figs ^ and ^1 
belong to different layers they move with different velocities. This is in sharp contrast with 
the behavior of sheared bulk suspensions where the layers always remain flat. 

In confined unsheared suspensions, the reservoir osmotic pressure, IT res , sets the volume 
fraction of particles in the gap. Then, the interplay between the particle osmotic pressure 
in the gap and the confined geometry determines the suspension structure J^, When 
the shear rate is large enough to produce non-equilibrium structures, the viscous stress, t v , 
dominates and determines the suspension structure. Therefore, the effective pressure arising 
from the shear stress [2^ must balance IT res to determine the volume fraction in the shear 
zone. If we start with a commensurate configuration and increase the gap keeping the layers 
flat, r v will decrease since there will be more fluid between the layers. To maintain the 
balance between t v and IL es the layers must increase their volume fraction. However, due 
to the constraints imposed by the confined geometry, sufficiently dense packing cannot be 
achieved with flat HCP layers. Instead, the system must adopt a buckled configuration. 
The shear stress forces the suspension to have translational invariance along the direction 
of particle motion. Thus the suspension cannot adopt a configuration which corresponds 
to an optimized 3-D packing. Instead, it must optimize its packing within the y'-z plane. 
Indeed, eliminating the fluid gaps and grouping the particles into triangles (Fig. EJ), we 
find the buckled structures resemble the densest 2-D packing of static hard disks under 
incommensurate confinement 




FIG. 2: Suspension structure in the y'-z plane. 
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To estimate t v , we consider a gap where the suspension forms flat HCP layers. Since 
the shear takes place in the lubricating fluid between the layers, the effective viscosity, 
rjeff, of the combined structure is: T] e ff = t)qD/ ^2 i l%, where k is the thickness of the ith 
lubricating fluid layer. The viscous stress is t v « / ~ffr] e ff. Intriguingly, the HCP layer 
structure introduces an ambiguity in the calculation since different flow configurations lead 
to different U. Nevertheless, we can calculate r v for two limiting flow configurations. In the 
straight flow configuration, particles move directly over the peaks of the HCP sheet below. In 
this case, we set k = when the peaks touch and the inter-layer distance is d. In the zigzag 
low configuration, particles move above the valleys formed by particles in the layer below 
^]. In this case, we set k = when the interlayer distance is dy/3/2. This is the minimum 
separation for flow configurations constrained to move without transverse displacements. 
By measuring the interlayer spacing for the different flows, we find that t T] « 6.0dyn/cm 2 
independent of the flow configuration. This stress must balance H res to ensure no net flux 
of particles between the shear zone and the reservoir. Simulations of unsheared disordered 
hard spheres [^J show n = l.ldyn/cm 2 when cf> = 0.61. Furthermore, in charged spheres 
n can easily reach six times this value j^J. Thus, the stress balance is consistent with our 
observations. 

This stress balance helps determine the non-equilibrium structure of the suspension for 
small D. We define D = D/d, set 7/ = 30s -1 , and, in Fig. El plot t v versus D for structures 
with up to 10 layers. The solid and dashed curves correspond to r v values calculated for 
straight and zigzag flow configurations respectively. The horizontal dash-dot line denotes 
Tires = 6.0 dyn/cm 2 . At each gap, the system must assume a structure where t v = TL res . 
Therefore, as D is reduced, a system initially in a straight flow configuration must increase 
the amount of zigzag with which the flat HCP layers move. For systems with fewer than 
eight layers, the stress curves for such configurations always reside between a solid curve 
at high D and a dashed curve at low D (Fig. EJ). However, as D is reduced further, the 
system must form a buckled structure with one fewer layer. Further reduction in D causes 
the amplitude of the buckles to decrease. Eventually the layers become flat and the system 
is described by the next solid stress curve with one fewer flat layers. 

We can further confirm this picture by accounting for the disappearance of the buckled 
state at large gaps. The D separation between the dashed curves in Fig. El is smaller than 
that of the solid curves. Consequently, the buckled state regions shrink with increasing D 
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FIG. 3: Transition mechanism predicting buckling in confined suspensions. The solid and dashed 
curves indicate t v versus the rescaled gap, D, for systems with up to 10 layers. The solid curves 
correspond to flat HCP sheets moving in the straight flow configuration. The dashed curves 
correspond to the zigzag flow configurations and delineate the border between the zigzag and 
buckled phases. With increasing gap, the buckled region becomes smaller and eventually disappears 
when the dashed and solid curves cross. The horizontal dash-dot line corresponds to the value of 
iT res . 

and eventually vanish when the curves cross. The curve crossing indicates that more than 
one flow configuration satisfies the pressure balance for gaps with more than 8 layers. To 
investigate this crossover in the experiments, we fix 7/ = 30s _1 and plot the separation 
between the top and bottom layers, AD t b, for the maximally buckled (open symbols) and 
straight (closed symbols) flow configurations with different numbers of layers (Fig. HJ. This 
measurement indicates the crossover appears when the system reaches 12 layers. The mis- 
match in crossover values suggests that in the straight flow configuration, the HCP sheets 
also move with some degree of registration so that the minimum interlayer separation is 
less than d. Indeed, even in the buckled state, where layers move with the least amount 
of zigzag, ping-pong-ball-models show the minimum interlayer separation is 0.94c?. This 
separation would lead to a crossover at 13 flat layers which is in excellent accord with our 
observations. 

An additional test of the model is obtained by observing how the slope of the linear fit to 
the buckled state data in Fig. 0] changes with shear rate. By varying 7/, the position where 
the dashed r n curves intersect the Il res line shifts. Therefore, the manner in which the D 
spacing between the t v curves varies predicts the shear rate dependence of the slope. For 
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FIG. 4: Separation between the top and bottom layers, AD^, for maximally buckled (open sym- 
bols) and straight (closed symbols) flow configurations with different numbers of layers. The 
measurement error is on the order of the symbol size. The shear rate 7/ = 30s -1 . A crossover of 
the curves occurs when the layer number is 12. 

0.3s -1 < 7/ < 30s -1 , the range of shear rates over which the measurements could be taken, 
the predicted slope is 0.95 at 7/ = 30s -1 and slowly approaches 0.87(~ Vo/2) as the shear 
rate is reduced. In the experiments, we observe a slope of 0.92 ± 0.05 at 7/ = 30s -1 which 
gradually decreases to 0.89 ±0.05 at 7/ = 0.3s -1 . Again, the measurements are in excellent 
agreement with the predictions. 

Our experiments clearly elucidate the complex interplay between the reservoir osmotic 
pressure and the shear stress in determining the unusual packings of confined colloidal sus- 
pensions subjected to high shear. Our results represent an important instance in which the 
non-equilibrium configurations of a sheared suspension can be determined. Similar effects 
will occur and must be considered in other flow geometries provided the shear stress dom- 
inates in the shear zone and the osmotic pressure dominates in the reservoir; for example 
such an interplay would be expected in a standard rheometer with a Couette geometry. This 
interplay may allow for tuning of the sheared suspension's rheological properties. Finally, 
since the observed structures are a consequence of sphere packing in confined geometries, 
such packing considerations may also impact the trends observed in the shear of very thin 
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